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^ ■ Aims. We investigate the non-Local Thermodynamic Equilibrium (non-LTE) line formation of neutral carbon in late-type stars 
' in order to remove some of the potential systematic errors in stellar abundance analyses employing C I features. 
f****) , Methods. The statistical equilibrium code MULTI was used on a grid of plane-parallel ID MARCS atmospheric models. 

Results. Within the parameter space explored, the high-excitation C I lines are stronger in non-LTE due to the combined effect 
of line source function drop and increased lower level population; the relative importance of the two effects depends on the 
particular combination of T e «, log g, [Fe/H] and [C/Fe] and on the analysed C I line. As a consequence, the non-LTE abundance 
Q ' corrections are negative and can be substantially so, for example ~ —0.4 dex in halo turn-off stars at [Fe/H]~ —3. The magnitude 
S—j \ of the non-LTE corrections is rather insensitive to whether inelastic H collisions are included or not. 

CZ) . Conclusions. Our results have implications on studies of nucleosynthetic processes and on Galactic chemical evolution models. 
When applying our calculated corrections to recent observational data, the upturn in [C/O] at low metallicity might still be 
present (thus apparently still necessitating contributions from massive Pop. Ill stars for the carbon production), but at a lower 
level and possibly with a rather shallow trend of ~ —0.2 dex/dex below [0/H]~ — 1. 
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1. Introduction 

Carbon is synthesised through the 3a nuclear reaction 
in stellar interiors (Burbidge et al. 1957; Arnett 1996; 
Pagel 1997). It is one of the most abundant elements in 
the Universe, together with hydrogen, helium, oxygen and 
neon. It can be found in a variety of compounds, so that 
it forms the basis of organic chemistry and of life (at least 
as we know it). Carbon is also important for dust forma- 
tion processes in the interstellar medium (e.g. Dopita & 
Sutherland 2003). Since it has a quite high ionization po- 
tential (11.26 eV), it remains mostly neutral in the photo- 
sphere of solar-type stars. Stellar spectroscopists thus rely 
on high-excitation permitted C I features and on the for- 
bidden [C i] line at 872.7 nm, as well as on various molecu- 
lar species (mainly CH and C2), to derive solar and stellar 
C abundances (e.g. Asplund et al. 2005a). 



A number of studies have investigated the carbon 
content of disk stars in the Galaxy (e.g. Andersson & 
Edvarsson 1994; Gustafsson et al. 1999), arguing that 
[C/Fe] 1 is slowly decreasing with time (i.e. with increasing 
metallicity) and probably mostly contributed by massive 
stars (MSs). In particular, Gustafsson et al. rule out a 
main origin of carbon in low-mass stars while they ar- 
gue that the contribution of intermediate-mass stars to 
the production of carbon in the Galaxy is unclear. Takeda 
& Honda (2005) have recently studied the CNO abun- 
dances of 160 disk F, G, and K dwarfs and subgiants with 
— 1 <[Fc/H]< +0.4 to look for any differences between 
the abundances of stars with and without planets. In par- 
ticular, when taking their calculated non-LTE corrections 
into account, they also found that [C/Fe] progressively de- 
creases with increasing metallicity for these galactic disk 
stars, by around ~ 0.2 dex per dex. A similar slope has 
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been found for [C/Fe] by Reddy et al. (2003, 2006) in 
their analysis of a number of elements in a large sample 
of F and G dwarfs in the thin- and thick-disk. In contrast 
with the above authors, Bensby & Feltzing (2006) find es- 
sentially no trend of [C/Fe] versus [Fe/H] for metal-rich 
thick- and thin-disk stars and conclude that the carbon 
enrichment at such metallicities is mainly due to low- and 
intermediate-mass stars (LIMSs), with massive stars still 
playing the role of main contributors at lower metallicities 
in the thick disk/halo. 

The halo stellar population is less well studied. It is 
still currently a matter of debate whether for halo stars 
carbon is produced primarily in LIMSs (e.g. Chiappini et 
al. 2003), or in massive stars (e.g. Akerman et al. 2004). 
It could well be that contributions from both of these 
sources operating on different timescales are important. 
Akerman et al. were the first to derive a previously unrec- 
ognized [C/O] upturn at low metallicities. However, they 
warned about possible non-LTE effects at low metallicity 
for the high-excitation C I lines they employed and cau- 
tioned that their results might be prone to uncertainties 
due to their assumption that C i-based abundances would 
be affected by the same non-LTE corrections as for Oi 
(which they estimated for their sample by interpolating 
between the values given in Nissen et al. 2002) and thus 
that the [C/O] ratio would be unchanged with respect 
to LTE. They mentioned the possibility that non-LTE ef- 
fects for the carbon and oxygen high-excitation lines might 
differ by an amount increasing with decreasing metallic- 
ity. At that stage, they were not able to exclude that the 
[C/O] values were being overestimated by ~ 0.2 dcx under 
the assumption of LTE. The idea that the upturn could 
thus be removed is in fact presented in some recent work 
(e.g. Takeda & Honda 2005; Fabbian et al. 2005; Asplund 
2005; Bensby & Feltzing 2006). Spite et al. (2005) study a 
sample of extremely metal-poor giants to investigate the 
abundances of C, N, O and Li. They find that in their 
"unmixed" stars (which should reflect the abundances in 
the early Galaxy) the [C /Fe] ratio is remarkably constant 
(at a level of ~ +0.2 — 0.3 dex) over the entire range 
-4.0 <[Fe/H]< -2.5 and that the [C/O] ratio is close to 
solar at the lowest metallicities. According to the studies 
of Akerman et al. (2004) and Spite et al. (2005), the up- 
turn in the [C/O] abundances of halo stars present at low 
metallicity could signal substantial early nucleosynthetic 
production of carbon in massive Pop. Ill stars. However 
it is worth noting that the inclusion of non-LTE/3D ef- 
fects plays a substantial role in testing such hypothesis. 
For example, 3D effects affecting the C and O abundance 
indicators used in the study by Spite et al. (2005) might 
reach ~ —0.6 dcx or more for CH, while they likely re- 
main around —0.2 dcx for the [O i] forbidden oxygen lines 
(Asplund 2005; see also Collet et al. 2006). This would 
then make their resulting [C/O] trend more or less flat 
at a level of ~ —0.6 dex down to very low [O/H] ratios, 
implying that resorting to high yields from massive zero- 
metallicity stars might be premature. 



On the theoretical side, low-metallicity fast-rotating 
models of massive stars (Meynet et al. 2006; see also 
Chiappini et al. 2006) predict important CNO yields 
through large mass loss during the red supergiant phase 
via stellar winds caused by efficient mixing and dredge- 
up enriching the stellar surface in heavy elements. These 
yields, especially if used in conjuction with a top-heavy 
IMF, would thus give large carbon abundances at low 
metallicity without introducing Pop. Ill stars. Carigi et 
al. (2005) have constructed various Galactic chemical evo- 
lution models with different C, N and O yields. The only 
models that are able to fit the carbon gradient are those 
including C yields that increase with metallicity due to 
stellar winds in massive stars and decrease with metallic- 
ity due to stellar winds in LIMSs. They argue that the 
fraction of carbon in the interstellar medium due to MSs 
and LIMSs would thus be dependent on time (with MSs 
dominating the C enrichment at early times, while around 
12 + log(0/H)~ 8 the LIMSs formed in the halo should 
be at the end of their evolution and they would start con- 
tributing with C ejecta in a comparable amount to that of 
MSs, see also Bensby & Feltzing 2006) and on galactocen- 
tric distance (LIMSs giving a larger contribution at large 
r Gc)) with similar contribution from both sources to the 
present carbon abundance in the interstellar medium of 
the solar vicinity. However, no attempt was made in that 
work to homogenize the sample of available observational 
data. Moreover, their "best" models are still not able to 
properly fit the N abundance values for metal rich stars 
and the C/O values at intermediate or low metallicity ob- 
tained from stellar observations. Gavilan et al. (2005) have 
recently suggested that a Galactic Chemical Evolution 
(GCE) model including yields from Woosley & Weaver 
(1995) for massive stars needs to be integrated with yields 
from LIMSs to properly reproduce the wealth of observa- 
tional data on the [C/Fe] trend down to [Fe/H]~ -2.4, 
suggesting that a downward revision of mass loss by stel- 
lar winds from massive stars is necessary and adding to 
the increasing evidence that LIMSs are significant con- 
tributors to the Inter Stellar Medium (ISM) carbon en- 
richment. Interestingly, Ballero et al. (2006) test several 
sets of yields for very massive zero-metallicity stars, using 
the two-infall chemical evolution model by Chiappini et 
al. (1997). They argue that the contribution of massive 
Population III stars to the Galactic evolution of carbon, 
nitrogen and iron would be negligible even at the lowest 
metallicities currently probed and that the existence of 
such objects could in any case not be proven or excluded 
on the basis of available data on known metal-poor halo 
stars, which likely formed as second-generation objects af- 
ter the explosion of primordial supernovae. 

Absorption systems give the opportunity to measure 
abundances in protogalactic structures at high redshift. 
For example, Becker et al. (2006) recently detected an 
overabundance of Oi absorption features in spectra of 
background high-redshift quasars (QSOs) towards the 
line-of-sight of Lyman-a forest clouds, deriving a mean 
< [C/O] >= -0.31 ±0.09. Moreover, Pettini (private com- 
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munication) and Erni et al. (2006) also find [C/O] 0.3 

in metal-poor Damped Lyman-a systems (DLAs) with 
[O/H]^ —2.5, in reasonable agreement with the values for 
Galactic stars (Akerman et al. 2004). These results are 
obviously precious in complementing Galactic studies at 
low [Fe/H]. In particular, Erni et al. argue that enrich- 
ment from a population of massive (10-50 M ) Pop. Ill 
stars exploding as core-collapse supernovae (SNell) and 
possibly as hyper-energetic (E> 10 51 erg) so-called hy- 
pernovae (HNe) is the most likely scenario, in contrast to 
models using non-zero metallicity progenitors or other ex- 
plosion mechanisms, such as super-massive (140-260 M ) 
pair-instability supernovae (PISNe) or Type la supernovae 
(SNela). Levshakov et al. (2006) derive a very low abun- 
dance of 13 C towards a quasar at z = 1.15. This does 
not support the enrichment of gas by fast-rotating, very 
low-Z massive stars (e.g. Meynet et al. 2006) and puts a 
tight bound to the possible amount of contributions from 
intermediate-mass AGB stars to the C nucleosynthesis in 
matter sampled by the line of sight to this system. 

The still open problem of the origin of carbon in the 
early Universe constitutes the main motivation for the 
present non-LTE study. The quality of stellar spectro- 
scopic observations possible at present, require us to make 
corresponding progress in models of stellar atmospheres 
and line formation processes in order to remove systematic 
errors. The impact on Galactic chemical evolution studies 
of such non-LTE calculations promises to be rewarding. 

2. C I non-LTE line formation 

2.1. Non-LTE preamble 

Spectral lines are powerful diagnostic tools that provide us 
with the capability of "remotely sensing" the physical con- 
ditions of the matter composing the particular atmosphere 
in which they form, via the analysis of the radiation we re- 
ceive from it. However, the problem of how these lines are 
actually formed when energy is transferred through the at- 
mosphere by radiation, is extremely complex (Mihalas & 
Mihalas 2000; Ruttcn 2003). The abundances of the chem- 
ical elements present in the atmospheres of stars are there- 
fore usually derived relying on the simplifying assumption 
of Local Thermodynamic Equilibrium (LTE) to derive the 
atomic level occupation numbers (Gray 1992). The visi- 
ble electromagnetic radiation we observe from stars comes 
from their photosphere, where radiation-matter interac- 
tion can play an important role. Radiation is in fact ab- 
sorbed and emitted by atoms and converted to different 
wavelengths via the interaction with the matter compos- 
ing the atmosphere; viceversa the atomic level populations 
are influenced by radiation. The LTE treatment assumes 
that this complex process is governed by a single temper- 
ature parameter. 

However, radiative rates often dominate over colli- 
sional ones in the atmospheric layers where the lines are 
formed and thus the LTE approximation can give mis- 
leading results. When LTE is not assumed a priori, the 



Saha-Boltzmann equation can no longer be used to deter- 
mine the level populations from local properties. A large 
number of radiative and collisional processes needs to be 
treated simultaneously to account for the coupling be- 
tween the populations of different atomic levels for the 
particular chemical element of interest. The major non- 
LTE effects in late-type stars have recently been reviewed 
by Asplund (2005). The main mechanisms are: 

— UV overionization: neutral minority stage elements 
(like e.g. Fei in the Sun) can be affected by overion- 
ization triggered by ultraviolet radiation (bound-free 
pumping). This is caused by a J;,-B„ radiation excess 
at short wavelengths (where J„ is the angle-averaged 
mean intensity and B u is the Planck function), with 
the radiation temperature exceeding the local electron 
temperature; 

— Resonance line scattering and photon losses: scattering 
from deep layers and photon losses can make the line 
source function drop below the local Planck function, 
resulting in non-LTE line strengthening; 

— Infrared overrecombination: occurs when J u drops be- 
low B„ and causes recombination to high-energy-lying 
states; 

— Bound-bound pumping: ultraviolet/optical radiation 
can cause pumping of bound-bound transitions, in- 
creasing the populations of upper levels; it can more- 
over cause overionisation (like in B I or Li i) , since other 
levels than the one involved in the transition can be 
directly influenced; and 

— Photon suction: caused by a recombination downflow 
from the first ionization stage, being for example the 
cause of the observed 12 /xm emission Mgi features in 
the Sun. 

2.2. Model atom 

Given the complex interaction of the different physical 
mechanisms in the line formation process, the accuracy 
of the atomic data needed as input is crucial to the final 
reliability of non-LTE calculations. 

Our carbon model atom contains 217 energy levels (207 
for C I, 9 for C II and one for C in) and is constructed tak- 
ing into detailed account a total of 650 radiative transi- 
tions (453 bound-bound and 197 bound- free). The neutral 
carbon levels are included up to an excitation potential 
of 11.13 eV (i.e. ~ 0.13 eV below the continuum, com- 
plete to a principal quantum number of n=9), the singly 
ionized levels up to 28.87 eV, while the Cm level is at 
35.64 eV. All radiative transition probabilities (f-values 
and photo-ionization cross-sections) are taken from the 
Opacity Project database TOPbase (Cunto et al. 1993; 
see also, specifically for carbon, Luo & Pradhan 1989 
and Hibbert et al. 1993). The necessary data for central 
wavelengths, radiative broadening and excitation poten- 
tial for the atomic energy levels were taken from the NIST 2 
(Martin et al. 1995; Wiese et al. 1996) database. 

2 http:/ /physics. nist.gov/PhysRefData/ 
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Our model atom also includes collisional excitations 
and ionizations by electrons and neutral hydrogen atoms. 
The cross-sections for forbidden transitions by electron 
collisions are derived with van Regemorter's (1962) for- 
mula (adopting an oscillator strength f= 0.01), while the 
impact approximation is employed for radiatively allowed 
transitions. One of the largest sources of uncertainty in 
current non-LTE studies is due to the unavailability of 
theoretical and laboratory collisional cross-sections for 
inelastic collisions with Hi. Hydrogen collisions in gen- 
eral act as a thermalizing source, bringing the balance 
of level populations closer to LTE. Especially in metal- 
poor stars, where neutral hydrogen atoms outnumber free 
electrons, such collisions should be of importance (e.g. 
Asplund 2005). We treated these collisions through the use 
of Drawin's formula (1968), as generalized by Steenbock & 
Holweger (1984) to the case of collisions between different 
particles. 

The Grotrian term diagram for our model atom is 
shown in Fig.^ When using this atomic model, we did 
not include fine structure splitting for the levels of interest 
here, which means that the level populations were com- 
puted with compound lines. We then expect each of these 
specific "grouped" lines to be stronger (and thus formed in 
higher atmospheric layers where the non-LTE effects are 
likely to be somewhat greater) than the corresponding ob- 
served splitted lines for the same term transitions. To be 
able to compare our results with observations, we have 
used a second model atom that includes fine splitting, re- 
distributing the converged population densities according 
to statistical weight. Table Ogives the relevant details for 
the lines included in this smaller atomic model (see also 
Fig.nj. Since one could be concerned that non-LTE effects 
may be overestimated using this approach 3 we have car- 
ried out tests which indicate that only relatively small dif- 
ferences are experienced in terms of abundance corrections 
when adopting the more consistent approach (the results 
remaining within ±0.05 dex of our standard results at all 
metallicities, the difference becoming even more negligible 
towards [Fe/H]~ —3). Since this marginal difference can 
be either positive and negative at the different metallici- 
ties, it is unlikely to introduce any systematic bias when 
using our computed non-LTE abundance corrections. In 
these tests, the radiative and collisional cross-sections in 
the splitted level case were reduced accordingly to the in- 
creased number of levels with respect to the unsplitted 
case, namely redistributing according to statistical weight, 
and a smaller 33-level atomic model was used to speed up 
the calculations, since we found that only small differences 
are obtained with respect to runs using the full 217-level 
atom (the results remaining within ±0.05 dex at all metal- 
licities). This is indeed a very interesting result and one 
which we plan to exploit in a future work on carbon non- 
LTE line formation using more advanced 3D atmsopheric 
simulation snapshots. The 3D non-LTE calculations in 

3 Because the fine-structure levels are not directly included 
in the whole problem 



Table 1. Details of the absorption lines included in the 
atomic model used after converged populations were ob- 
tained. A total of 11 transitions were treated in detail. 
Wavelength of absorption line, upper and lower level iden- 
tification, excitation potential for lower level and oscillator 
strength are respectively given for each of the transitions 
included 
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fact become feasible in a reasonable amount of time with 
this relatively simple atom. 

The standard 217-level atom we use in this study is 
sufficiently complex to also include transitions commonly 
used as solar diagnostics. For example, our study confirms 
that the 538.0 nm line experiences small non-LTE effects 
(< 0.05 dex in a solar-like model, as also found by Asplund 
et al. 2005a): solar observations could be used to test the 
non-LTE modelling. This feature becomes too weak be- 
low [Fe/H]~ —1 and thus is not useful for determining C 
abundances at low metallicity. In this work we have in- 
stead paid particular attention to the high-excitation per- 
mitted Cl lines around 910.0 nm, which are detectable 
down to low metallicity and have been used by Akerman 
et al. (2004) . Our aim is to study departures from LTE for 
these lines, with a focus on the resulting C abundances in 
metal-poor halo stars once the non LTE corrections are 
taken into account. 

2.3. Method 

The statistical equilibrium code MULTI 4 (Carlsson 1986, 
1992) was used to simultaneously solve the radiative trans- 
fer and rate equations for the C I non-LTE spectral-line 
formation. Complete redistribution over Voigt line profile 
has been assumed. A grid of plane-parallel ID MARCS 
(Gustafsson et al. 1975; Asplund et al. 1997 and subse- 
quent updates) model atmospheres with varying T e g, log g 
and [Fe/H] has been employed in the calculations. The 
atmospheric parameter grid covers the range of late-type 
stars, comprising a total of 168 models with 4500 < T e s 
< 7000 K, 2.0 < log g < 5.0 [cgs] and -3 < [Fe/H] < 0. 
In all cases, unless noted otherwise, the microturbulence 
velocity has been assumed to be £ = 1.00 km/s. An inher- 

4 Available at |http://www.astro.uio.no/~matsc/mul22/| 
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Fig. 1. Left panel: Grotrian term diagram for our 217-level (207 for C I, 9 for C n and one for C in) carbon model. For 
clarity, only the 453 bound-bound radiative transitions for Cl are shown here. Right panel: Grotrian term diagram 
for the smaller carbon model atom used to compute the non-LTE abundance corrections for the Cl lines around 
830.0 — 960.0 nm of particular interest here. The model includes 10 Cl levels and 1 Cn level. The bound-bound 
transitions in the triplet system generate from splitted levels in this model, producing three and six lines around 
963.0 and 910.0 nm respectively 



ent assumption in our study is that the departures from 
LTE for carbon do not influence the given atmospheric 
structure. This is reasonable in the sense that carbon is 
not an important direct opacity nor a main electron con- 
tributor, so that the atmospheric opacity (and thus the 
equilibrium structure of the atmospheric model) will not 
change significantly even when large non-LTE corrections 
are found. 

The standard solar C abundance adopted in our non 
LTE radiative transfer calculations has been logec = 
8.40, very close to the recently updated value (Asplund 
et al. 2005a, b). For each model atmosphere, the carbon 
abundance was varied to cover the range —0.60 < [C/Fe] 
< +0.60 in steps of 0.2 dex. We have repeated all calcu- 
lations for three different choices of the scaling factor Sh 
for inelastic collisions (Sh = 0, 10 -3 and 1). 

3. Results 

3.1. Specific stars: Sun, Procyon, HD140283 and 
G64-12 

Initial tests were carried out for a few models represent- 
ing specific stars (Sun, Procyon, HD 140283, G64-12), to 
check the effect of changes in different atomic properties 
and identify the main driving non-LTE effects for different 
stellar parameters. 

The high-excitation C I lines studied here are much 
stronger in Procyon (T cff = 6350 K, log g= 3.96, [Fe/H] = 



0.00, Allcnde Prieto et al. 2002) than in the Sun 
(5780/4.44/0.00), so that the higher line formation height 
in Procyon results in larger non LTE effects. In the Sun, 
non-LTE abundance corrections vary between ~ —0.10 
and ~ —0.25 dex for the lines of interest (when neglect- 
ing H collisions), while they rise to between ~ —0.25 and 
<~ —0.50 dex in Procyon. In both stars the drop of the 
line source function below the Planck function drives the 
departure from the LTE approximation. However, since 
the level populations deviate more from the LTE case 
when moving outwards in the atmosphere, in the case of 
Procyon the non-LTE abundance corrections are more se- 
vere. 

The low-metallicity of HD 140283 (5690/3.67/-2.40) 
implies that the lines are weak and thus formed relatively 
deep in the atmosphere. Naively one would then expect 
the non-LTE effects to be quite small. Still, the non-LTE 
abundance corrections in this star vary between ~ —0.20 
and ~ —0.35 dex, i.e. larger than in the Sun. This is 
mainly due to the lower log g and thus densities imply- 
ing less efficient collisions. The lower [Fe/H] also results 
in smaller electron pressures at a given optical depth and 
therefore again less important collisions. The level popu- 
lations start to deviate relatively deep in the atmosphere 
of HD 140283, causing a line source function drop and an 
increased line opacity (because of the lower level of the 
transition getting overpopulated) , both effects acting at 
the same time to strengthen the line. The non-LTE line 



6 



D. Fabbian et al.: Cl non-LTE spectral-line formation in late-type stars 



formation in G64-12 (6511/4.39/-3.20/) is similar to HD 
140283: the lower log g and [Fe/H] and higher T e g com- 
pared to the Sun all work in the same sense, thus implying 
even more pronounced non-LTE abundance corrections, 
between —0.2 and —0.45 dex, mainly due to an overpop- 
ulated lower atomic level (opacity effect making the line 
stronger than in LTE). 

We have investigated the sensitivity of the non-LTE 
calculations to the adopted electron and hydrogen colli- 
sional cross-sections. We ran calculations with a few dif- 
ferent values for the parameter Sh regulating the H col- 
lisions efficiency. We found small abundance variations of 
the order of ~ 0.1 dex when changing Sh even by large 
factors, suggesting that the resulting non-LTE corrections 
are not particularly sensitive to the exact choice. This also 
guided us in what values to adopt (Sh = 0, 10~ 3 and 1) 
when successively running full sets of calculations for the 
whole grid of stellar models described below. We checked 
the effect of electron collisions by artificially decreasing 
their cross-sections by a factor of 10. Also in this case, the 
effects were small and non-negligible only at low metallic- 
ity where they reach ~ ±0.03 dex at most. 

We finally tested the influence on the non-LTE correc- 
tions when increasing the microturbulence velocity from 
1.0 to 1.5 km/s. A higher value for this parameter, intro- 
duced in ID atmospheric models as an effective (though 
not perfect) substitute for convective motions, is expected 
to affect mostly the stronger saturated lines, since the 
atoms will absorb photons at slightly shifted wavelengths 
(where the flux is higher) compared to the line center. 
The stronger lines will thus be strengthened and lower 
abundances should be obtained from them. By appropri- 
ately adjusting the microturbulence parameter, it is pos- 
sible to remove trends of derived abundances with exci- 
tation potential or line strength which appear when ve- 
locity fields are ignored. As expected, the lines' equiva- 
lent widths (both LTE and non-LTE) increase slightly, 
by only a few mA when increasing the microturbulence to 
£ = 1.5 km/s. However, the resulting non-LTE abundance 
corrections experience insignificant variations, decreasing 
by < 0.01 dex. 

3.2. Non-LTE effects 

The departure coefficient (3 for an atomic level is defined 
as (3 = N/Nitfe, where N and Nlte are the atomic pop- 
ulations in non-LTE and LTE respectively, for the level 
considered. The ratio of departure coefficients for different 
atomic levels enters the line source function 5 Sj,, so that 
(when neglecting stimulated emission, as should be possi- 
ble for the lines considered here) identical values of (3 U and 
Pi (for the lower and upper level of a transition), result 

5 More explicitly, assuming complete redistribution and in 
the Wien regime, the ratio of the line source function to the 
Planck function can be expressed in terms of the ratio of depar- 
ture coefficients for the upper and lower levels of a transition, 
as S l v /B u » p u /Pi 



in a Planckian line source function. This means that even 
when both levels in a transition have populations that de- 
part considerably from the Boltzmann expectation, it is 
possible that the net effect is a line source function close 
to LTE, as long as the departures are the same for the 
two levels. However, the effect of line opacity must also be 
considered when the lower level population increases or 
decreases with respect to the LTE approximation, since it 
will affect the depth of line formation. 

The leftmost panels in Fig.|3show, for a representative 
case (T e ff= 6000 K, log g= 4.50 and collisional efficiency 
set to Sh = 10 -3 ), the typical departure coefficients f3 for 
the upper and lower levels of the C I line at ~ 910.0 nm. 
The figure highlights how at solar metallicity the domi- 
nant non-LTE effect is the line source function dropping 
below the Planck function, mostly due to the decreased 
population of the upper level compared to LTE, while at 
low metallicity the line source function is close to LTE, 
in particular at the much deeper formation layers in such 
stars, i.e. closer to the continuum optical depth; the lo- 
cation of T continuum = 1 and T Une cen ter = 1 is marked 
in the figure. The non-LTE effects are found to be large 
even for the most metal-poor stars in our grid. To explain 
this, we used the contribution functions defined by Magain 
(1986) to evaluate the contribution of the different atmo- 
spheric layers to the formation of the spectral-line depres- 
sion. Especially at low metallicity we found contributions 
from layers much higher up than the canonical line optical 
depth unity to be important. The non-LTE line strength- 
ening effect can then be explained in such case, in terms 
of increased line opacity: the lower level of the transition 
gets overpopulated with respect to LTE in those layers, as 
seen in the departure coefficients plot for this metallicity 
(Fig.EJl, causing the still large non-LTE corrections found 
in this case. In stars with [Fe/H]~ —1, the two effects just 
described (line source function drop and line opacity in- 
crease) work in unison, resulting in large non-LTE effects 
(see Fig.|HJ). 

Figs.|21and0]show results from "multi-MULTI" runs. In 
such calculations, all 650 radiative transitions and a subset 
of all collisional transitions (all collisional rates between 
the lowest 26 atomic levels) were perturbed in turn by a 
factor of 2 to study the impact on the C I lines of main 
interest here. The runs were carried out for a star with 
T e s= 6000 K and log g= 4.50 and metallicities from solar 
to [Fe/H]= —3. Since for every transition that is changed 
in the atomic model a full MULTI run has to be performed, 
the multi-MULTI runs are very time-consuming. We refer 
to Carlsson et al. (1994) for a detailed description of how 
such runs are carried out. 

Fig-IHl confirms that for the ~ 910.0 nm C I feature, 
the driving non-LTE effect at solar metallicity is the line 
source function drop. The ~ 166.0 nm resonance line from 
the ground state (2p 2 3P - 3s 3 P°) is an extremely opti- 
cally thick transition and cannot drain the population of 
its upper level (which is the lower level of the ~ 910.0 nm 
transition), which will thus stay close to LTE (/?/ ~ 1) 
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Fig. 2. Departure coefficients from LTE (left panels), line source and Planck function comparison (middle panels) and 
resulting line profile (right panels) for the Cl line at ~ 910.0 nm shown for different [Fe/H] for runs with our larger 
atomic model. The atmospheric models all have T c g— 6000 K, log g— 4.50 [cgs] and [C/Fe]= 0, with [Fc/H] going from 
0.0 (top) to —3.0 (bottom). The hydrogen collision efficiency has been set to Sh = 10~ 3 . For the departure coefficients 
the solid curve represents the lower level and the dotted one the upper level for the transition of interest. In the 
middle panels the line source function (solid line) and Planck function (dashed line) are normalized to the continuum 
intensity. Location of T cont i n uum = 1 (solid line at logio T500 ~ 0) and Tu ne center = 1 (hi LTE and non-LTE, dashed 
and solid line respectively) are shown by vertical lines in the left and middle panels. In all cases, the non-LTE line 
profile (solid line) is stronger than in LTE (dashed line) 



where the line is formed due to sufficiently strong cou- 
pling to the ground state through electron collisions (see 
Fig.|3J). The upper transition level instead gets underpop- 
ulated in non-LTE (f3 u = TV/JVlte < 1), mainly through 
photon losses in the line itself, see Fig.03 Since the ratio 
Pu/Pi determines the line source function, this explains 
why the latter drops below the Planck function (S l v < B u ) 
around solar metallicity (Fig.|2J) , causing the line strength- 
ening effect. The effect of collisions (mainly in the triplet 
system to which the line studied belongs) is instead to 
weaken the line, tending to restore the LTE populations, 
however the line source function drop still dominates, driv- 
ing the non-LTE line strengthening. At this metallicity, 



the ~ 910.0 nm Cl line will therefore be most sensitive 
to processes that affect the upper level population, thus 
changing the line source function. There is very little or 
no coupling with the singlet system at this metallicity, so 
that only a small influence from collisions in the singlet 
system is seen. Photoionization is not important either. 



For metal-poor stars, the line source function is close to 
LTE in the deep layers where the line is formed. Indeed, 
even though both levels tend to get overpopulated with 
respect to LTE, the ratio of their departure coefficients 
stays close to the LTE value {fi u /fli ~ 1) where the line 
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Fig. 3. Results from multi-MULTI runs for a star with T c g= 6000 K and log g— 4.50 and solar metallicity are shown 
for the Cl line around 910.0 nm. The effect on the line strength of multiplying each of the different rates for radiative 
transitions in turn by a factor of 2 is shown in the left panel. The right panel shows the effect when varying the 
collisional rates individually by the same amount. Solid lines indicate positive contribution (strengthening), dashed 
lines mean negative contribution (weakening) from that particular transition to the strength of the line studied. The 
thickness indicates the relative change in the non- LTE equivalent width (W pcrt -W rcf )/W rof between the perturbed 
and reference case, as in the figure legend, while the arrows indicate the net flow (evaluated at Tune center = 1-0) for 
the transitions 



is formed 6 . However, the onset of significant departures 
from LTE in the level population is now shifted to deeper 
atmospheric layers than for the solar case. Moreover, the 
overpopulation of the lower level in these layers makes 
the line form further out in the atmosphere (where devi- 
ations in the level populations are more dramatic) than 
in LTE, thus they are stronger than in that case. This 
opacity effect causes the line strengthening to be sen- 



It is important to remember that, as mentioned, these weak 
lines form with a significant contribution from a tail of atmo- 
spheric layers located much higher up than one would expect 
in the simple Eddington-Barbier approach, as we discovered 
by looking at their contribution functions (Magain 1986). For 
example, at T e «= 6000 K, log <?4.50, [Fe/H]= —3, the average 
height of formation actually moves out to logio T500 ~ —0.7, 
with a significant contribution to the absorption profile still 
coming from layers as high up as logio T500 ~ — 2. Also note 
the steep increase of the level population at low metallicity, as 
one moves outward in the atmospheric model (Fig.|5J 



sitive to processes affecting the lower level population. 
Even though electron collisional coupling to the ground 
state still tends to restore LTE population in the lower- 
level of the ~ 910.0 nm transition also at low metallic- 
ity, the latter level now gets severely overpopulated across 
the atmosphere, due to the fact that the total sum of the 
net radiative rates dominates with respect to collisional 
ones for this transition (see Fig.[!3) and that it cannot de- 
populate efficiently and will act as a meta-stable level in 
detailed balance with the ground state. Moreover, in the 
low-metallicity case, the line strength is influenced by pro- 
cesses in the singlet system much more than at solar metal- 
licity, i.e. the coupling between the triplet (to which the 
lines studied belong) and singlet systems becomes impor- 
tant (Fig. 2} and starts to influence the line strength more 
and in a more complex way, both through radiative rates 
and collisional transitions. So, while at solar metallicity 
the magnitude of the non-LTE effects is mainly governed 
by collisions and radiative transitions between levels in the 
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triplet system itself, with decreasing metallicity both the 
intersystem (triplet-singlet) collisions and the processes in 
the singlet system (collisions as well as radiation) gain im- 
portance with respect to processes in the triplet system. 
Electron collisions, especially those between levels with 
similar energy, activate the coupling with processes in the 
singlet system and thus for example play an important 
role in determining the population balance of the lower 
level of the transition of interest. The use of accurate col- 
lisional rates is therefore important. Even though at low 
metallicity the lines studied are weak, thus formed deeper 
in than for solar metallicity (Fig.[2J), and one could naively 
assume the non-LTE effects to be small, it is the relative 
strengthening of the line with respect to LTE that justi- 
fies the still large non-LTE corrections we found for low- 
metallicity stars. The main non-LTE effect as one moves 
to low metallicity is thus an increased line opacity arising 
from an overpopulated lower level, with the contribution 
coming from the line source function effect still present but 
becoming less important with respect to the solar metal- 
licity, because both the lower and upper levels of the tran- 
sition are now similarly overpopulated in the line-forming 
layers. 

Note that in Figs .[31 and 0] the upward or downward 
direction of the rates (represented by arrows) is evalu- 



ated at Tu ne cen ter = 1-0. In reality such rates can have 
a complicated behaviour with atmospheric depth. For ex- 
ample, for the left panel (solar metallicity, radiative case) 
in Fig.[3 the rate for the transition from level 2s 2p 3 3D° 
at ~ 7.95 eV to 2s 2 2p 3p 3P 2 at - 8.85 eV (i.e. from level 
6 to 14, in the notation adopted in Fig.^) is upward, but it 
is found to become downward just one depth-point below 

7~Une center 

= 1.0. The net result of this is that the up- 
per level population density for the C I transition around 
910.0 nm actually decreases when increasing this f- value, 
which explains why such line becomes stronger (a line 
source function effect, because (3 U /Pi becomes less than 
unity, thus making S;, <B y ). 

In Fig.[S]we plot the net radiative and collisional rates 
for the transitions that most influence the lower and up- 
per level of the transition at ~ 910.0 nm. One for example 
sees that at solar metallicity the population balance of the 
upper level (3p 3 P2, labelled as level 14 in the figure, as 
per notation adopted here of increasing level number with 
increasing energy, with level being the ground state), 
shown in the middle left plot, is set by radiative losses to 
the lower level of the transition itself (3s 3 P°) balanced by 
collisions from this same level and from level 3p 3 Pi which 
has energy close to that of the upper level. The line will 
thus get stronger if any of the levels feeding into the up- 
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Fig. 5. The net rates for radiative (no markers) and collisional (with markers) transitions are shown in the upper left 
plot for the lower level of the transition at ~ 910.0 nm in a model with T e g= 6000 K, log g= 4.50 and solar metallicity. 
Only the rates for the transitions important for the level balance are shown, marked with different symbols. Their 
sum is also marked with crosses and confirms that the rest of the levels play a minor role here. The horizontal line at 
dn/dt= marks the expectation that the total sum of all rates is equal to zero. The levels in the legends are numbered 
accordingly to our atomic model (see Fig.^), i.e. increasing with excitation potential. The middle left plot similarly 
shows the rates for the upper level of the transition, while the bottom left plot shows the total sum of the net radiative 
and collisional rates respectively (positive=upwards). The plots on the right show the same quantities, but for a model 
with [Fe/H]= —3. Note the different scale for the vertical axes in the plots with different metallicity 



per level get underpopulated, conversely it will get weaker 
if we increase the population of any of these. This pic- 
ture confirms what is seen in Fig.|21and described above, 
namely that collisions generally have a weakening effect 
on the line, by tending to restore LTE populations. The 
sum of all net radiative rates is predominant in magni- 
tude with respect to that of collisional rates, only slightly 
in the solar-metallicity case, while it strongly dominates at 
low metallicity (bottom panels of Fig.^J). The effect is to 
drive the balance away from LTE, with a decreased upper 
level population in the solar metallicity case mainly due 
to photon loss in the line itself (the lower level population 



remaining close to LTE due to collisional transitions out 
of that level which tend to mitigate non-LTE effects by 
contrasting the effect of radiative transtions from higher 
levels) and an increased population for both levels at low 
metallicity (the population flow from higher levels into the 
lower level of the transition causing its overpopulation and 
thus driving the departure from LTE). As line transitions 
turn out to be important, in different ways, both at so- 
lar and low metallicity, our choice of including as many 
levels and accurate radiative rates in the atomic model 
employed (see Sect. l2.2j) . making it more complete than 
those used in some previous studies by other authors (e.g. 
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Fig. 6. Non-LTE abundance corrections over the grid of temperatures and gravities explored, with varying metallicity 
(as indicated in the insets), for the neutral carbon absorption line at 909.4834 nm (used e.g. by Akerman et al. 2004 
as one of the abundance indicators for their determination of C in halo stars). The corrections refer to the case 
with [C/FeJuE = +0.4, typical of metal-poor halo stars and have been obtained from calculations with fine splitting 
included in the atomic model. The four leftmost panels show the corrections we obtained when neglecting hydrogen 
collisions, while the rightmost panels represent the results with Sh = 1. Large non-LTE corrections are present at 
high temperature and low gravity, in the various cases 



Stiirenburg & Holweger 1990; Rentzsch-Holm 1996a, b), 
is particularly well suited. 

Our results underline that the non-LTE effects for the 
carbon lines considered have different causes in different 
stars. It is thus paramount to carefully evaluate such ef- 
fects on a star-to-star basis, in stellar abundance studies. 

3.3. Non-LTE abundance corrections 

The C I non-LTE line strength for the parameter grid 
explored is generally larger than in LTE. The non-LTE 
abundance corrections reach down to as much as —0.8 dex 
and more for particular combinations of atmospheric pa- 
rameters (notably at T cS > 6000 K, log g= 2.00, [Fe/H]= 
— 1.00, with corrections of up to ~ —1.0 dex). For each set 
of atmospheric parameters, the resulting non-LTE abun- 
dance corrections become more severe for stronger lines, 
similarly to what is found for the solar case (Asplund et al. 
2005a) . This can be explained intuitively because stronger 
lines will be formed higher up in the atmosphere where de- 
partures from LTE in the level populations are often more 
significant. 

Several studies (e.g. Kaulakys 1985; Fleck et al. 1991; 
Bclyaev et al. 1999; Belyaev & Barklem 2003) now seem to 
indicate that Drawin's classical recipe (1968) might over- 
estimate the efficiency of H collisions by about three or- 
ders of magnitude. In view of the corresponding uncer- 
tainty of our final results due to this treatment of the 



excitation and ionization by H collisions, we have chosen 
to perform full sets of calculations for the whole parame- 
ter grid, applying Drawin's formula with different choices 
of the scaling factor Sh for separate runs, to test how 
much this affects our derived non-LTE abundance correc- 
tions. We have adopted Sh=0, 10 -3 and 1 respectively 
and have found a relatively small sensitivity of the re- 
sulting non-LTE corrections to the value of this parame- 
ter. The non-LTE abundance corrections obtained for all 
stellar models in our grid, from calculations adopting the 
two extreme choices of Sh are shown in Fig.[B]for the car- 
bon absorption line at 909.4834 nm. The plot refers to 
results obtained when adopting [C/Fe]= +0.4, which is 
the typical LTE value for the most metal-poor halo stars 
in the Akerman et al. (2004) sample. As expected, we find 
slightly larger non-LTE abundance corrections when set- 
ting a low efficiency for H collisions. In all cases, stars 
in the higher temperature, lower gravity range appear to 
be particularly affected, with very large non-LTE effects. 
The corrections peak around a metallicity [Fe/H]~ —1. As 
seen in the figure, when neglecting inelastic H collisions or 
assuming Drawin's recipe (1968) for their treatment, the 
non-LTE corrections only vary marginally (< 0.1 dex dif- 
ference between corrections obtained in the two cases). 
Since our choices of Sh should bracket the true values, 
it is reassuring that the corresponding uncertainty in the 
non-LTE corrections is small. 
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Fig. 7. The non-LTE abundance corrections are plot- 
ted versus line strength for the C I lines around 
963.0 and 910.0 nm. The models have T cS = 6000 K , 
log g= 4.0 and [Fe/H]= -3, -2, -1 and (for which dif- 
ferent symbols are used in the figure: filled squares, dia- 
monds, triangles and circles respectively). Here [C/Fe]= 
and Sh = 10~ 3 have been assumed. For each absorption 
line, the non-LTE corrections at different metallicities are 
connected to underline the effect of the metallicity change 



At the lowest metallicity, the non-LTE effects are still 
important, with Alogec ~ —0.35 dex for typical halo 
stars parameters (T e s ~ 6000 K, log g= 4.00) and a choice 
of Sh = f -3 . When we instead adopt Drawin's recipe as 
is (parameter regulating collisions with Hi set to Sh=1), 
the resulting non-LTE corrections remain important and 
are only slightly smaller (still amounting to ~ —0.25 dex 
at the lowest metallicity for typical halo stars parame- 
ters) than the previous case. Fig. shows the large non- 
LTE corrections obtained for a representative set of atmo- 
spheric parameters (T c ff=6000 K and log g=4.00). They 
indeed reach down to ~ —0.45 dex at [Fe/H]= —1 and 
are still very substantial (A log ec ~ —0.35 dex) even for 
weak lines at [Fe/H]= —3. Note that in these plots, only 
corrections obtained with [C/Fc]= +0.00 are shown. For 
larger [C/Fe] the line formation is shifted outwards and 
hence the non-LTE corrections become even more severe. 

Our analysis aims at removing (some of the) poten- 
tial systematic errors in stellar abundance analyses. With 
this in mind, we provide an IDL routine that interpolates 
our predicted non-LTE abundance corrections to give the 
result for arbitrary stellar parameters and carbon abun- 
dances within our grid. 

4. Comparison with previous non-LTE studies 

Only a small number of previous studies have investigated 
the non-LTE formation of carbon lines in late-type stars 
(Sturenburg & Holweger 1990; Takeda 1994; Paunzen et 

7 Available via FTP from: 

ftp:/ /ftp. mso.anu.edu.au /pub /damian /corr_nlte/c / 



al. 1999; Takeda & Honda 2005; Asplund et al. 2005a). 
With the exception of Takeda & Honda (2005), these have 
all been restricted to a few specific stars rather than cov- 
ering a large parameter space as we did in our study. 

To determine the carbon abundances of their sam- 
ple of disk stars, Takeda & Honda (2005) use the C I 
505.2/538.0 nm permitted lines. These abundance indica- 
tors were also used by Asplund et al. (2005a) for a deter- 
mination of the solar photosphcric carbon abundance. In 
both cases, very small non-LTE effects were found for such 
lines. Wc confirm that these features have small negative 
non-LTE abundance corrections (A log ec no larger than 
~ —0.05 dex for solar-type stars). To reanalyse the metal- 
poor halo stars available in the literature, Takeda & Honda 
carried out non-LTE calculations down to [Fe/H]= —4 in 
their work, to determine the departure coefficients for the 
levels relevant to the same high-excitation features we also 
study in detail. They used a carbon model atom with 129 
Cl levels (Takeda 1992) for an extensive parameter grid 
similar to ours but adopting [C/Fe]= and ±0.3 dex only 
in the calculations. For the collisional rates due to neutral 
hydrogen atoms, they also adopted Drawin's formula in 
the generalized version by Steenbock & Holweger (1984), 
with the default choice Sh = 1 for all calculations. We 
find a generally good agreement between the two non-LTE 
studies. However for a few models with high T e g and low 
log g and [Fe/H] our non-LTE abundance corrections, al- 
though very large (Alogec — —0.8 dex), are not as severe 
as claimed by Takeda & Honda (see our Fig.|HJl. 

We suspect that the reason for this difference is an 
extrapolation error in the work by Takeda & Honda, pos- 
sibly due to their use of only [C/Fe]= and ±0.3. For the 
halo stars available in the literature, they also find rela- 
tively large (|Alogec| £ 0.3 dex), negative, metallicity- 
dependent non-LTE corrections below [Fe/H]= —2, those 
values however being still ~ 0.05 dex smaller than what 
we obtain at such low metallicities. 

Sturenburg & Holweger (1990) found that C non- 
LTE abundance corrections in the Sun are typically small 
(;S —0.05 dex on average when adopting a scaling factor of 
1/3 to Drawin's formula for the H collisions). For the lines 
we also study here, they explained the non-LTE effects in 
the solar photosphere mainly in terms of lower level over- 
population, causing the line absorption coefficient to ex- 
ceed the LTE value. We underline that the driving effect at 
solar metallicity is instead that of the line source function 
dropping below the Planck function, i.e. both the lower 
and upper levels of the transition need to be considered. 
The upper level underpopulation is actually more signif- 
icant than the lower level (slight) overpopulation in the 
line-forming region, driving the line source function drop 
in Sun- like stars (see our Fig.|21). In any case, a detailed 
comparison between the results for the high-excitation C I 
lines of interest shows a very good agreement, the non- 
LTE abundance corrections amounting to ~ —0.1 dex for 
those features in both studies. 

Using the same model atom as Sturenburg & Holweger, 
Tomkin et al. (1992) find comparably small non-LTE ef- 
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Fig. 8. The non-LTE abundance corrections we found in 
this work (vertical axis scale) for the 909.4834 nm absorp- 
tion line in a range of atmospheric models are compared 
to the corresponding results in Takeda & Honda (2005), 
for the case Sh = 1. The solid line marks the one to one 
relation, while the dotted lines include a region where the 
difference between the corrections in the two studies is 
limited to < 0.15 dex 



fects for the Cl ~ 910.0 nm lines (Aloge c < -0.15 dex) 
also in metal-poor halo stars. However, they adopted very 
large efficiencies for the H collisions, in order to reproduce 
the observed line strengths using the high solar carbon 
abundance expected at the time. As explained, evidence 
is now accumulating that such high values for H collisions 
might be unrealistic by several orders of magnitude, so 
that the non-LTE corrections for the high-excitation C I 
features of interest should be expected to be in reality 
more significant than in that study. 

Using an extended version of the atomic model by 
Sturenburg & Holweger, Rentzsch-Holm (1996a, b) in- 
vestigated carbon non-LTE abundance corrections in A- 
type and related stars. A detailed comparison with that 
work is unfortunately not possible, since the two calcu- 
lations only overlap for our highest temperature value 
(7000 K). However, for that common temperature value, 
a good agreement is found when comparing the resulting 
non-LTE corrections for the high-excitation C I lines; at 
log g= 3.5, Alogec reaches ~ —0.45 dex in both stud- 
ies. For a given stellar model, we too find increasing non- 
LTE corrections with line strength, see the trend indi- 
cated by identical symbols in Fig.[7| this effect has al- 
ready been explained by Sturenburg & Holweger (1990) 
as due to the numerical enhancement, in the flat part of 
the curve-of-growth, of the abundance correction derived 
from equivalent width changes due to non-LTE, compared 
to the abundance correction derived in the linear part of 



the curve-of-growth. When decreasing gravity or increas- 
ing effective temperature (see our Fig. El note that for 
T e s> 5500 K, the non-LTE abundance corrections tend 
to become more insensitive to temperature changes when 
fixing metallicity and gravity, except at very low metallic- 
ity), generally larger non-LTE corrections result for both 
studies. Another similarity between these two works - and 
also with that of Sturenburg & Holweger 1990 - is that 
the four lowest levels of Cl (not involved in the high- 
excitation transitions of interest here) in our model are 
found for almost all cases in the parameter space we study 
to have Boltzmann populations throughout the model at- 
mosphere, due to strong collisional coupling between the 
neutral ground state and the next three lowest-lying lev- 
els. This ensures that the forbidden [Cl] line at 872.7 nm 
(transition between the two lowest-lying levels in the sin- 
glet system) is in general not affected by departures from 
LTE. The levels involved in this transition only start to ex- 
perience non-LTE effects at [Fe/H]= -3, for T cS > 6500 K 
and higher (and increasing with decreasing gravity too). 
The two levels tend to experience rather similar depar- 
tures from LTE. Rentzsch-Holm (1996b) similarly found 
that increasing T g for the higher temperature and metal- 
licities range in that study resulted in undcrpopulation of 
these levels: they showed this can transmit to other levels 
through line transitions and lead to global underpopula- 
tion of Cl around ~ 12000 K. In our case, both the high 
temperature and low metallicity act to increase the magni- 
tude and the atmospheric depth of the onset of departures 
from LTE, because the UV flux can more effectively un- 
derpopulate the atomic levels involved in the transition. 

5. Galactic chemical evolution of carbon 

Equipped with our non-LTE calculations, we have re- 
visited the observational study of C and O in halo stars by 
Akerman et al. (2004). Those authors found that [C/O] 
seems to rise for [0/H]< —1; which could suggest high 
carbon yields from zero-metallicity, so called Population 
III stars. They however assumed that the non-LTE abun- 
dance corrections for the two elements would be the same. 
We have calculated specific non-LTE corrections for their 
sample of 34 dwarf halo stars. Since in their analysis those 
authors used the same atmospheric models as we employ 
here, our derived non LTE corrections can be directly ap- 
plied to their LTE abundance results. Considering their 
estimate for random errors in equivalent widths measure- 
ments and stellar parameters determination - which in 
turn affects their resulting LTE abundances - and the un- 
certainty in our non-LTE calculations (~ 0.1 dex), we con- 
sider it feasible to study the C and O trends. Moreover, 
the [C/O] ratio is essentially insensitive to changes in the 
stellar parameters, given that the C and O high-excitation 
abundance indicators used are affected in a similar way. 
We are thus confident that our results are a step forward 
with respect to previous works. 
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Table 2. The IDs, atmospheric parameters (effective temperature, gravity and iron content) and the LTE abundances 
for oxygen and carbon for the sample of low-metallicity stars from Akerman et al. (2004) are shown in column 1, 2, 
3, 4, 5 and 7 respectively. The corresponding LTE [C/O] values are in column 10. Our adopted non-LTE abundance 
corrections for oxygen and carbon (obtained when neglecting collisions with Hi atoms) and the resulting [O/H] and 
[C/O] non-LTE abundance ratios from this study are shown in column 6, 8, 9 and 11 respectively. All values from the 
Akerman et al. study have been corrected, where appropriate, for our different choice of the carbon and oxygen solar 
abundances (logec Q = 8.39 and logeo = 8.66) 
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*Note that the O I non-LTE abundance corrections for the five most metal-poor stars have been revised with respect to those 
given in Akerman et al. (2004), as explained in the text 

* The C and O abundances of G016-013 have been revised by +0.4 dex with respect to those in Akerman et al. (2004), due to 
the fact that the wrong (hotter) stellar model for HD160617 was accidentally used in that analysis instead of the correct one for 
this star (P.E. Nissen, private communication). This affects [O/H] (and [C/H]) but leaves [C/O] intact for this star. We trust 
that the stellar parameters for this star in Akerman et al. (also given here) have been correctly derived in their analysis 



Table 2 lists the results obtained for the C and O abun- 
dances. Fig. [5] shows the trend with atmospheric param- 
eters of the carbon non-LTE corrections we obtained for 
the sample in Akerman et al. (2004). A clear increase in 
the non-LTE abundance corrections (becoming approx- 
imately 0.3 dex more negative at [Fe/H]~ —3 than at 
[Fe/H]~ —1) is seen as one moves to low metallicity, which 
is mostly due to the higher effective temperatures of the 



sample stars towards the metal-poor regime (T c g values 
increase by approximately 500 K or more for stars below 
[Fe/H]~ —2, compared to the higher-metallicity stars in 
the sample). This is still in agreement with the fact that 
in general, when only decreasing the metallicity of the 
atmospheric model to [Fe/H]~ —3, the non-LTE abun- 
dance corrections are roughly of the same order than at 
solar metallicity (see Figs. El and 0). A similar increase be- 
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Fig. 9. The non-LTE abundance corrections (empty triangles: Sh = 0; filled triangles: Sh — 1) for the 34 halo stars 
in the sample by Akerman et al. (2004) as a function of T c s (a), log g (b) and [Fc/H] (c) 
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Fig. 10. Trends of [C/Fe] vs. [Fe/H] (left panel) and [O/Fe] vs. [Fe/H] (right panel) in Milky Way halo (squares) and 
disk stars (circles, Bensby & Feltzing 2006, from [C i] and [O i] lines free from non-LTE effects). The LTE values from 
Akerman et al. (2004), corrected for our different choice of carbon and oxygen solar abundances, are represented as 
empty squares in the two panels, while the non-LTE abundances are shown as filled squares. The non-LTE results 
are those obtained when neglecting collisions with hydrogen. Note that the [O/Fe] non-LTE values in the right plot 
for the five most metal-poor stars already include the new O I non-LTE abundance corrections we present here (see 
text) 



low [Fc/H]^ —1 in the non-LTE abundance corrections 
of Akerman et al.'s sample of halo stars was found by 
Takeda & Honda (2005), however when we adopt their 
same choice for the H collision efficiency (Sh = 1) and mi- 
croturbulence (£ = 1 km/s), we still find corrections that 
are more negative than theirs by ~ 0.05 — 0.10 dex for the 
lowest-metallicity stars in that sample. 

Previous studies of Milky Way disk stars (e.g. 
Gustafsson et al. 1999) found a slowly decreasing [C/Fe] 
with increasing metallicity (slope of —0.17 dex per 
dex). However, very recently, Bensby & Feltzing (2006), 
analysing a sample of thin and thick disk stars and us- 
ing [Cl] and [Oi] lines immune from non-LTE effects, 
found a more or less flat trend from around solar metal- 
licity down to [Fe/H]~ —0.9. We find a roughly solar 
[C/Fe] (or slightly higher when adopting Sh = 1) extend- 



ing much further down to [Fe/H]~ —3.2 (Fig. llO|) . albeit 
with a rather large and - at least in part - likely real scat- 
ter. The plot shows how the LTE [C/Fe] values in the 
Akerman et al. study have a large scatter around a value of 
~ +0.35 dex, with a possible hint for a slight increase with 
decreasing [Fe/H]. This trend is removed when correcting 
for our increasingly large (as one moves down in metal- 
licity) non-LTE corrections for the stars in their sample. 
In Fig.^]we also show [O/Fe] as a function of [Fe/H] for 
the halo stars of Akerman et al. As for the [C/Fe] ver- 
sus [Fe/H] plot, the disk stars of Bensby & Feltzing are 
included too. The original results of Akerman et al. have 
here been corrected for our revised O I non-LTE correc- 
tions at the lowest [Fe/H] as well as for our different so- 
lar abundances (Asplund et al. 2005b). The [O/Fe] trend 
in metal-poor stars is still a hotly debated topic and a 



16 



D. Fabbian et al.: Cl non-LTE spectral-line formation in late-type stars 



0.2 

0.0 
-0.2 



o 



H. -0.4 



-0.6 



-0.8 



■1.0 - 



* ■ 

LP815-43 



□ 
□ 
□ 



-3.0 -2.5 -2.0 



□ □ 
_ ■ u 

"m a □ 
■ ■ ■ 



■ J ■ 

□ □ 



1.5 -1.0 -0.5 0.0 0.5 
[0/H] non _ LTE 



0.2 - 
0.0 - 
-0.2 - 
-0.4 - 
-0.6 - 
-0.8 - 





• 






□ 


■ M » 


T - 





.o h 

-3.0 



□ LP815-43 □ 



□ 
□ 
□ 
4 



■ J 3 
■ □ r 

£ D J D 

H ' 

e i 

□ 



□ □ • ■ 

□ 



□ 



-2.5 



-2.0 



1 .5 

ro/Hi 



-1 .0 



0.0 



0.5 



Fig. 11. Trend of [C/O] vs. [O/H] in Milky Way halo and disk stars. The upper panel represents the new observational 
trends obtained when our non-LTE calculations are applied to the data in Akerman et al. 2004 (using our different 
choice of carbon and oxygen solar abundances and assuming negligible collisions with hydrogen). The [O/H] values 
already take into account our more negative non-LTE corrections for oxygen at low metallicity. The empty and filled 
squares represent the [C/O] LTE and non-LTE values respectively. The lower panel uses the same notation, with 
results now calculated accounting for hydrogen collisions for carbon (with a choice of Sh = 1)- In both panels, the 
starred symbol indicates the [C/O] value we believe is more realistic for LP8 15-43 (see discussion in text) 



full account of this problem is beyond the scope of the 
present paper; we refer the reader to recent discussions 
such as Nissen et al. (2002) and Asplund (2005) and ref- 
erences therein. We note however that Fig. suggests a 
roughly flat trend at a level [0/Fe]~ +0.6 with no in- 
dication of a linear trend with slope of ~ —0.3 dex/dex 
as claimed by Israelian et al. (1998) and Boesgaard et al. 
(1999). Asplund & Garcia Perez (2001) have argued based 



on 3D hydrodynamical model atmospheres of metal-poor 
stars that O abundances derived from OH lines as used 
by Israelian et al. and Boesgaard et al. can be severely 
overestimated in ID (see also Collet et al. 2006). 

A few peculiar cases in the two plots in Fig.llQI should 
be noted: 

- The C and O LTE abundances for G16-13 ([Fe/H]= 
—0.76) are both increased by 0.4 dex with respect to 
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Akerman et al. (2004), due to the fact that a too 
hot stellar model was used in their analysis for this 
star (P.E. Nissen, private communication). This moves 
the star from the lower to the upper extreme of the 
distribution in the [C/Fe] and [O/Fe] versus [Fe/H] 
plots. Because the scatter in both plots is rather large, 
this does not significantly affect the overall resulting 
trends. Moreover, since both the C and O abundances 
are increased by the same amount, the star maintains 
the same [C/O] as in Akerman et al. (demonstrat- 
ing how insensitive [C/O] is to the model used), while 
[O/H] increases by 0.4 dex: the overall [C/O] trend is 
unchanged too, again because of masking by the rel- 
atively large scatter at the higher-metallicity extreme 
of the distribution, to which the star is shifted 
— The two most extreme cases in Fig.^5] are LP815- 
43 (with a low [0/Fe] no „-LTE = 0.37) and G20- 
8 (showing very high [C/Fe] non _LTE = 0.31 and 
[0/Fe] non _ LTE = 1.02). We note that for the O I triplet 
equivalent widths in LP815-43, Nissen et al. (2002) 
give different values than those in Akerman et al. 
(2004). The smaller line-to- line scatter in the former 
determination seems to better correspond to the vari- 
ation in gf- values; moreover, the S/N for this star is 
higher in Nissen et al.'s spectrum. These authors derive 
log eo = 6.71 compared to log eo = 6.61 in Akerman 
et al. (2004). It is thus possible that the oxygen abun- 
dance derived in the latter work needs to be increased 
by up to 0.1 dex, with the errors in the equivalent 
width determinations likely attributable to e.g. con- 
tinuum placement problems. The case of G20-8, which 
shows anomalously large [C/Fe] and [O/Fe] values (by 
~ 0.3 dex), seems to indicate that its temperature 
might be underestimated by as much as 400 K, possi- 
bly due to uncertainties in the reddening 

The O I non-LTE abundance corrections are taken 
from Akerman et al. (2004), with the exception of the 
five most metal-poor stars. Akerman et al. interpolated 
the non-LTE corrections of Nissen et al. (2002) as a func- 
tion of O I equivalent width, leading them to deduce a 
non-LTE correction of —0.10 dex for those five stars. We 
have carried out non-LTE calculations in the same man- 
ner as described by Nissen et al. Our corrections for the 
five most metal-poor stars are more negative than found 
by Akerman et al. (2004), varying between —0.18 and 
—0.30 dex. Moreover, while - with respect to [C/O] - we 
find that the non-LTE corrections for carbon and oxy- 
gen almost balance each other at higher metallicity, at 
[0/H]< —1 they become more negative for carbon. While 
the Oi non-LTE corrections are typically ~ —0.2 dex, 
and —0.3 dex at most at [Fe/H] ~ —3, the correspond- 
ing Cl non-LTE corrections are ~ 0.15 dex lower still 
at these low [Fe/H]. We adopt slightly different solar C 
and O abundances than done by Akerman et al. 8 , which 
also slightly modifies the resulting [C/H], [O/H] and [C/O] 

8 Namely, logec = 8.39 and logeo = 8.66 (Asplund et al. 
2004; Asplund et al. 2005b), as compared to logec = 8.41 



values. Had we used the solar values determined with a 
ID MARCS model atmosphere instead, [O/H] would have 
been lower by -0.05 dex and [C/H] higher by +0.04 dex 
only, for all stars. Clearly, such a small and constant off- 
set is not in any way crucial in determining the general 
trends of [C/Fe] and [O/Fe], because all stars would shift 
by the same amount. More importantly, the [C/O] values 
remain essentialy unchanged. Indeed, we recover a similar 
trend to that in the work by Akerman et al. for the [C/O] 
abundance ratio at low metallicity, even though at a lower 
level (see Fig. HT|) . We note that the similarity is due to 
the fact that the carbon non-LTE corrections in our study 
are more negative at low metallicity than those we esti- 
mate for oxygen, by approximately the same amount for 
all of the metal-poor stars below [O/H]^ —1.5. This gives 
a more or less constant downward shift of those [C/O] 
values as compared to the Akerman et al. study, who had 
assumed for carbon that the same non-LTE corrections 
would applicable as for oxygen. The revised (more nega- 
tive) O I non-LTE abundance corrections we use for the 
five most metal-poor stars in the sample partially com- 
pensate the large non-LTE corrections we find for Cl at 
such low metallicities. These increased non-LTE correc- 
tions for oxygen also shift those stars to correspondingly 
lower [0/H] non -LTE values with respect to the original 
study by Akerman et al., "stretching" the low- metallicity 
part of the plot with respect to their work. 

In summary, the adoption of the recent values of solar 
abundance and the application of our carbon and oxy- 
gen non-LTE corrections has the net effect of shifting 
down the [C/O] values of the low-metallicity stars by 
similar amounts, preserving the possible [C/O] upturn 
at [0/H]< -1, although at a lower [C/O] level. We do 
urge that this rising trend should not be over-interpreted 
at this stage. First of all, if our suspicion is confirmed 
that the oxygen abundance has been underestimated by 
~ 0.1 dex by Akerman et al. for LP815-43 - which is the 
star with the highest [C/O] ratio at low metallicity - it 
would imply that this star will accordingly have a lower 
[C/O] , bringing it closer to the rest of the low-metallicity 
points and giving in the end a less significant indication for 

the claimed upturn. The [C/O] trend below [O/H] 1 

might thus have a gradual increase with a shallow slope 
of ~ —0.2 dex/dex. Moreover, in view of the large depar- 
tures from LTE we found for Cl in ID, it is important 
to carry out 3D non-LTE calculations. Furthermore, our 
tests have revealed that the O I non-LTE effects at these 
very low [Fe/H] may be more uncertain than previously 
thought due to large sensitivity to poorly known collisional 
cross-sections. We are indeed planning to investigate this 
in an upcoming work. Finally, the T c g scale for halo stars 
is still under debate (e.g. Mclendcz et al. 2006). The claim 
of carbon production from high-mass Pop. Ill stars is thus 
still uncertain in our opinion. 



and logeo = 8.74 in Akerman et al. (2004). We believe the 
values we adopt are best estimates since they are averaged over 
different lines and they account for 3D/non-LTE effects 
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Spite ct al. (2005) used different abundance indicators 
(namely, the G band of CH and the forbidden [O i] line at 
630.0 nm) to extend the study by Akerman et al. down to 
[Fe/H] ~ —4 ([O/H]^ —3), again deriving an increasing 
[C/O] ratio reaching a near solar level (<~ —0.2 dcx) to- 
ward lower metallicity. However, they did not account for 
3D effects which should affect the carbon molecular lines 
they employed more than the [O i] line they used. Even 
though those authors studied extremely metal-poor halo 
giants, for which 3D effects are not well known, one can 
use the results for turn-off stars at similarly low metal- 
licity as an indication, so that the corresponding values 
given in Spite ct al. might turn out to be overestimated 
by ~ 0.3 dex or more (Asplund 2005) and the [C/O] trend 
could still be roughly constant or only slowly increasing. 

6. Conclusions 

We have carried out non-LTE spectral-line formation 
calculations for C I absorption features. The compu- 
tations have been performed on a grid of 168 ID 
MARCS stellar model atmospheres with varying param- 
eters (4500 <T cff < 7000 K, 2 <log g< 5, -3 <[Fc/H]< 0) 
and -0.60 <[C/Fe]< +0.60). The main results are: 

— We generally find large, negative non-LTE abundance 
corrections for high-excitation neutral carbon permit- 
ted lines, i.e. LTE abundances are too large. At so- 
lar metallicity, the dominant non-LTE effect is that of 
the line source function dropping below the Planckian 
mostly due to under-populated upper atomic levels in 
the transitions, but at low metallicity it is that of line 
opacity due to overpopulated lower levels 

— We find substantial non-LTE abundance corrections 
for metal-poor halo turn-off stars: typically ~ —0.4 dex 
at [Fe/H]^ —3; these results are relatively immune to 
the adopted H collisional efficiency 

— Using our resulting non-LTE corrections, we revisit 
the observational study of Akerman ct al. (2004) and 
find that the [C/O] trend at low metallicity might still 
be suggesting an upturn indicating possible production 
of carbon through nucleosynthesis in Pop. Ill stars. We 
however warn that more observational data is needed 
- together with detailed non-LTE calculations for oxy- 
gen, to understand the magnitude of the abundance 
corrections applicable for this element at metallicities 
[Fe/H]< —2.5 - in order to reduce the still large uncer- 
tainties and finally settle this issue 

Without doubt one of the necessary steps for obtain- 
ing reliable stellar abundances in future large scale sur- 
veys will be to perform a vast quantity of detailed non- 
LTE studies for an increased number of different chem- 
ical elements and stellar parameters. Very many precise 
atomic input data for collisional and radiative processes 
have to be known for realistic results. The use of sophis- 
ticated techniques for quantum-mechanical computations 
have led to accurate atomic radiative data fortunately 



becoming available in recent years to the astrophysical 
community, thanks to such large scale efforts as e.g. the 
Opacity and IRON Projects. We underline the necessity 
for more theoretical and laboratory studies aimed at ob- 
taining such accurate atomic data and in particular at 
improving our knowledge on the correct treatment of the 
still poorly understood collisions with electrons and with 
hydrogen atoms. The effect of using 3D hydrodynamical 
model atmospheres should be investigated as the 3D non- 
LTE effects could prove to be even stronger than in ID, 
due to the presence of large temperature inhomogeneities 
and to the different temperature structure in the optically 
thin atmospheric layers with respect to ID hydrostatic 
models, especially for the low-metallicity regime (Asplund 
et al. 1999). 
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